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Abstract 

We consider the interacting holographic dark energy with new in- 
frared cutoff (involving Hubble parameter and its derivative) in non- 
flat universe. In this context, we obtain the equation of state pa- 
rameter which evolutes the universe from vacuum dark energy region 
towards quintessence region for particular values of constant parame- 
ters. It is found that this model always remains unstable against small 
perturbations. Further, we establish the correspondence of this model 
having quintessential behavior with quintessence, tachyon, K-essence 
and dilaton scalar field models. The dynamics of scalar fields and 
potentials indicate accelerated expansion of the universe which is con- 
sistent with the current observations. Finally, we discuss the validity 
of the generalized second law of thermodynamics in this scenario. 
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1 Introduction 



Dark energy (DE) is one of the most attractive and active fields in modern 
cosmology due to the indications of accelerated expansion of the universe 
through type la Supernovae [1]. Observational data like CMBR [2], large 
scale structure [31 H], gravitational lensing surveys |5] and galaxy redshift 
surveys |6] also favor this phenomenon. However, the identity of DE is still 
ambiguous and various models have been suggested to know its nature. Its 
simplest candidate is the cosmological constant but it suffers two well-known 
problems, i.e., "fine tuning problem" and "cosmic coincidence problem" [7]. 

In order to understand the nature of DE phenomenon, various dynamical 
DE models have been proposed which can be characterized by the equation 
of state (EoS) parameter u. The holographic dark energy (HDE) is one of 
the emergent dynamical DE model proposed in the context of fundamental 
principle of quantum gravity, so called holographic principle [8]. It is de- 
rived with the help of entropy-area relation of thermodynamics of black hole 
horizons in general relativity which is also known as the Bekenstein-Hawking 
entropy bound, i.e., S ~ M^L'^, where 5* is the maximum entropy of the sys- 
tem of length L and Mp = {8tvG)~^ is the reduced Planck mass. Using this 
relation, Cohen et al. [9] argued that the vacuum energy (or the quantum 
zero-point energy) of a system with size L should always remain less than 
the mass of a black hole with the same size due to the formation of black 
hole in quantum field theory. Hsu [TU] and Li [TT] formulated this statement 
in mathematical form as 

which is known as HDE density, m is constant, L is the infrared (IR) cutoff 
and G is the gravitational constant. 

The density of HDE describes the connectivity between ultraviolet (UV) 
and IR cutoffs which represent the bounds of energy density and size of 
the universe, respectively. The HDE model suffers the choice of IR cutoff 
problem. Li [TT] proved that Hubble as well as particle horizons are not 
compatible with the present status of the universe while the future event 
horizon is the best candidate for non-interacting HDE with suitable param- 
eter m. It is argued fl2] that HDE with future event horizon plagued with 
causality problem (why should we calculate the current value of DE density 
with the help of future event horizon of the universe?). This problem moti- 
vated people [I2]-[I1] to modify the IR cutoff as a function of the Ricci scalar 
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or generalized form of the Ricci scalar. Further, observational analysis has 
also been done for these types of HDE models [TJ]. 

The stability against small perturbations is a well-known procedure to 
check the viability of a DE model. For this purpose, the sign of the square 
of the speed of sound, = plays a key role - its negativity represents 
instability and vice versa [7]. It was shown [TB] that Chaplygin and tachyon 
Chaplygin gases are stable as > 0. However, for the holographic [T7], age- 
graphic p!8j and QCD ghost DE [19] models, < 0, i.e., they are classically 
unstable. 

The scalar-field dark energy model then can be considered as an effective 
description of this holographic theory. The reconstruction of HDE in terms 
of scalar fields has been discussed widely. The scalar fields (which naturally 
occur in particle physics such as string theory [2D]) are used as a possible can- 
didate of DE. In this scenario, a large number of models have been proposed 
including quintessence, phantom, K-essence, tachyon, ghost condensates and 
dilatonic DE [20] [21] etc. Granda and Oliveros [22] formulated the scalar 
field models for HDE by using new IR cutoff [13] (called new HDE (NHDE)) 
in flat FRW universe. Karami and Fehri [23] generalized this work to the 
non-fiat universe. Recently, Sheykhi [23] has constructed the quintessence, 
tachyon, K-essence and dilatonic DE scalar field models for interacting HDE 
with Hubble horizon as an IR cutoff in fiat universe. The generalized sec- 
ond law of thermodynamics (GSLT) was also discussed for HDE model with 
different IR cutoff for non-fiat universe [25j. 

It is believed that an early inflation era provides a fiat universe, but 
this consequence is only true if the number of e-foldings is very large [26] . 
The data of first year WMAP analysis favors the non-fiat scenario of the uni- 
verse [27] . Also, different observational data provided the evidence about the 
contribution of spatial curvature to the total energy density of the universe 
[28]-[33]. It was shown that the parameterizations of the dark energy mod- 
els admit the non-flat universe by implying compatible observational data 
[35] . Recently, Lu et al. [36] have used type la supernovae, baryon acoustic 
oscillations, CMBR and observational Hubble data and obtained the value 
— 0.00131q;qq4q for fractional energy density due to curvature. It would be 
interesting to study the universe with a spatial curvature. 

In view of above discussion, we extend the work of Sheykhi [23] by taking 
the generalized form of IR cutoff of HDE in non-flat universe. We also check 
the validity of the GSLT in this scenario. The paper is organized as follows: 
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Section 2 contains the discussion of the evolution and instabihty of interact- 
ing NHDE in non-flat universe. Section 3 is devoted for the reconstruction 
of scalar field models of interacting NHDE while section 4 investigates the 
validity of GSLT. The last section provides the summary of our results. 



2 New Holographic Dark Energy 

In this section, we manipulate the expressions of EoS parameter and speed 
of sound for NHDE interacting with dark matter (DM) in non-fiat FRW 
universe 

ds" = ~df + a\t)[ ^ + r^de^ + sme^d(l)% (1) 

1 KV 

Here a{t) is the cosmic scale factor which measures the expansion of the 
universe and k = —1,0,1 represents the spatial curvature indicating the 
open, fiat and closed universes, respectively. The corresponding equations of 
motion are 

^' + 4 = i(P™ + PA), (2) 

Qi O 

H + H^ = -liPm + pA + ^PA), (3) 

o 

here we assume SnG = 1 as well as dust like DM. Also, pm and pa are DM 
and DE densities, p\ is pressure due to DE, H denotes the Hubble parameter 
and dot represents derivative with respect to time. We can rewrite Eq.([3]) in 
the form of fractional energy densities as 

fi^ + = 1 + ^k, (4) 

where 

O _ o _ Pa _ k 



3^2' - 3^2' - ^2^2 

We assume the NHDE density in the following form [13] 

PA = 3{pH^ + XH), (5) 

where p and A are positive constants. 

The interaction of DM and NHDE leads to the equation of continuity in 
the form of two non-conserving equations as 

Pm + SHpm = T, (6) 
Pa + 3H{pa+Pa) = -T. (7) 
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Figure 1: Plot of wa versus x for NHDE. 



We choose the following form of the interaction term T 

T = ^^Hprn, (8) 

here d"^ is the interacting constant. Using Eqs.([6]) and ([8]), we obtain DM 
density as follows 

where x = Ina and is an integration constant. Inserting Eqs.(l5]) and ([9]) 
in Eq.(l2]), we obtain the differential equation 

' 

where E'^ = (0) represents the present value of the parameter and b is 
an integration constant. This has the solution 

/i-A-1 2(^1 -1)- 3(1 -d^)X ^ ■ 

The initial condition = 1 yields a; = 0, E = 1, and hence we have 

5=1 ^ + ^^"^0 (12) 

//-A-1 2(/i-l)-3(l-rf2)A- ^ ' 
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Figure 2: Plot of versus x for NHDE. 



With the help of Eqs.© 
of interacting NHDE as 



and (fTTl) . we obtain the evolution parameter 



2d^n 



mo 



_3{fi-X-l] 
(3A-2/x + 2)e- 



2(// - 1) - 3(1 - d2)A^ ~3A 



-2(^-1)3: 



(/i-A-1) 



X (2(/i-l)-3(l-t/2)A)-ia„„e 



-3{l-d^)x 



+ be' 



(13) 



We plot the EoS parameter u\ versus x with respect to different well-known 
choices of interacting parameter (P, i.e., 0, 0.04, 0.08 shown in Figure 1. Also, 
we assume the NHDE parameters as /u = 1.198, A = 0.195 and the current 
values of Qko = 0.01, I^Ao = 0.73. We see that the EoS parameter translates 
the universe from vacuum DE region towards quintessence region. 
The speed of sound is given by |7j 



p 
P 



p 



(14) 



where prime means differentiation with respect to x. Using Eqs.(l5]), ([7]), f lTT]) 
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and fll4p . it follows that 



X 



X 



X 



2(/i - x)nko ^_2^ 



6rf2(i - d^)a 



mo 



_3(/i - A - 1 
(3A - 2/i + 2 

V 



2(^- 1) -3(1 -(i2)A 



-3(l-(i2)a; 



+ 26(// - 1) 



-2()j-l)3: 
A 



-2(^- A)fi 



L (/"-A-l) 



/CO g-2x 



+ (2/i - 3(1 - rf=^)A) 



3(1 - d2)(2(/. - 1) - 3(1 - d')X)-'n^,e-'('-''> - ^^^L_A 



(15) 



which provides the speed of sound for interacting NHDE. This is shown in 
Figure 2 for the same constant parameters as given above. In this scenario, 
the speed of sound remains negative with the increase of interacting param- 
eter. This shows that the NHDE is unstable just like non-interacting HDE 
model with future event horizon |17] . 



3 Reconstruction of New Holographic Scalar 
Field Models 

Here, we provide the correspondence of the interacting NHDE with quintessence, 
tachyon, K-essence and dilaton field models in non-flat universe. 

3.1 New Holographic Quintessence Model 

The ordinary scalar field is governed by quintessence which is minimally 
coupled with gravity. The energy density and pressure of the quintessence 
scalar field are defined as [2T] 



P, = l'P' + Vi<l>), p^ = ^<j>'-V{<P), (16) 

where 0^ and V{(fi) are termed as kinetic energy and scalar potential, respec- 
tively. The EoS parameter for this model becomes 

^ 2V-(0) 
+ 2^(0) ■ 
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Figure 3: Plot of (j) versus x for quintessence model. 



For the correspondence between NHDE and quintessence scalar field, we set 
Pq = Pa and Pg = p\. Consequently, Eq. (fT6!) yields 



6^ = 3H^ 



2{fi - x)nk, _ 

3(/i-A-l) 



2(^- 1) -3(1 -d2)A 



26(/i — 1) -2(^-1)^ 

e A 



3A 



(17) 



3(/i-A-l) 



2(2(^-l)-3(l-rf2)A) 



6(3A — u + 1) -2(m-i)^ 
+ 3A ' ' 



(18) 



The kinetic energy can also be written as 



2ip-m^^.2. _ {2p-3{l-d')X-2d')nmp ^.,a_,2), 



3(^ - A - 1 



26(/i — 1) -2(m-1)3: 
P A 

3A 



2(^- 1) -3(1 -d2)A 



-2x 



[(/i-A-i: 



(2// - 3(1 - rf')A) 



(2(// - 1) - 3(1 - d')A)-^l]„„e-=^(^-'^')" + 6e- 



-2()j-l)3: 



(19) 



This equation cannot be solved analytically to obtain (/) due to its com- 
plicated nature. To get insights, we solve it numerically and plot against x 



8 




Figure 4: Plot of V versus x for quintessence model. 
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by choosing initial condition, 0(0) = 0, while keeping the remaining param- 
eters same as in the previous section. Figure 3 shows that the scalar field 
increases (and hence the kinetic energy cj? of the potential decreases) with 
the passage of time. The potential V versus x and are shown in Figures 
4 and 5, respectively, indicating the decreasing behavior. The quintessence 
potential in terms of shows large value at the present epoch which repre- 
sents accelerated expansion of the universe. This is consistent with the result 
shown through phase space analysis [21] that the exponential potential for 
the scalar field contains the attractor solutions describing the accelerated 
expansion of the universe. However, in the later time, the universe remains 
in the accelerated phase because quintessence potential goes to positive and 
non-zero minima but kinetic energy goes to zero. 

3.2 New Holographic Tachyon Model 

The tachyon model, originated from the string theory, has been suggested 
to explain DE scenario. It has an interesting feature that a rolling tachyon 
interpolates the EoS parameter between —1 to 0. Also, the tachyon model 
is the best candidate for infiation at high energy. Many attempts have been 
made to formulate reliable cosmological models with the help of different self- 
interacting potentials [37]. However, the effective Lagrangian for this model 
is defined as pi] 

L = -y(0)^i + a^0a^0, 

where V^(0) represents the tachyon potential. This scalar field has the fol- 
lowing energy and pressure 

Pi = ^ffi=, p, = -F(0)^/rr^ (20) 
y^l-02 

and the EoS parameter is 

= 0' - 1. (21) 
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The correspondence between NHDE and tachyon model is obtained for pt 
Pa and pt = Pa which lead to 



_3(/i - A - 1) 
2b(n — 1) -2(^i-i)x 



{2p-3{l-d^)X-2d^)a 



2(2(/i- 1) -3(1 -d2)A) 
^''~^^^'°-e-'^-{2p-3il-d')X) 



lip-X-l] 



X {2{p - 1) - 3(1 - (i')A)-^fi^„e-3(^-'^')" + 6e A 



2()j-l)3: 



(22) 



y(0) 



3H', 



{2fi-3{l-d^)X)a 



+ be' 



[{p-X-1) 

2(tJ.-l)x 



2{p- 

ii^-x)nk, ^_2^ 

3(/i - A - 1) 



mo g-3(l-d2)a 



l)-3{l-d^)X 
-2d' Q„,,e-'^'~'''>{2{p-l) 



- 3(l-d^)A)-i + A(3A_2// + 2)e"^ 
3A 

From Eqs.l ITT]) and ( l22l) . we obtain the kinetic energy term as 



(23) 



+ 



1 

2b{p 



2Jji-^X)9^^_2^ 
3(^ - A - 1) 



(2// - 3(1 - rf2)A - 2(i2)^] 



rm_^-Z{l-(f)x 



3X 



1) -2(m-1)x 

— e ^ 



2(^- 1) -3(1 -d2)A 

(2/x-3(l -d2)A) 



(/i - A)fifco.g-2r. 



_(/i- A-1 



X (2(^-1)- 3(1 -d2)A)-il]™oe 
X (/i- A-l)-ie-2^ 



+ 6e' 



-2(^-1)3: 
A 



ko 



2fi^„e-=^(i-'^')^ , -2(M-i). 



2(/i - 1) - 3(1 - c/2)A 



.(24) 



The plot of the tachyon field (p versus x is given in Figure 6 which shows 
that the scalar field (j) increases with the passage of time and becomes steeper 
for increasing the interacting parameter d'^. However, the kinetic energy of 
the tachyon potential decreases and approaches to zero in the future. With 
this behavior, Eq.f l2ip indicates the vacuum evolution of the universe. The 
tachyon potential shows oscillation initially about its maxima (it increases 
with the increment of d^) but approaches to zero in the later time as shown 
in Figure 7. For the later time, its rapid decrease from maxima gives inverse 
proportionality to the scalar field. This type of behavior corresponds to 
scaling solutions in the brane-world cosmology 



11 





Figure 7: Plot of y versus x for tachyon model. 

3.3 New Holographic K-essence Model 

The K-essence model is different from quintessence scalar field model in the 
sense that it evolutes the universe in the accelerated expansion era. It is 
originated from the idea of K-inflation which was used to describe the infla- 
tion of the early universe at high energies |20] . This model has been used as 
an alternative candidate of DE which yields interesting results of scaling and 
attractor solutions [521110] ■ This model is described with a scalar field having 
non-canonical kinetic energy. The generalized form of scalar field action is 

m 
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where p{(f), x) shows the pressure density as a function of potential and 
X — \ i*^- The corresponding energy density and pressure are 

Pfc = y((/>)(-X + 3x'), Pk^vm-x + x^). (25) 



where V{(f)) represents the scalar potential of K-essence model. The corre- 
sponding EoS parameter is 



1-X 
l-3x' 



(26) 



Equating = Pa and Pk = Pa for the correspondence between NHDE and 
K-essence model, we obtain 



X = 



4(^-A)Qfe„ ^ _,, _ (2/x-3(l-rf^)A + 2rf^)l^^„ ^_3(,_,.), 

3(^ - A - 1) ^ 



b , , , -2(tJ.-l)x 

+ — (6A - 2// + 2)e^^ 
3A 



2(^-1) -3(1- d2)A 

e - - {2n - 3(1 - cZ^)A 



- X)Qko -2x /o., Q^l 



(^-A-1) 



+ 6rf2)(2(^ - 1) - 3(1 - c/2)A)-il^„„e-3(i-'^')- -h -(4A -2/^ + 2) 

A 



-2(n-l)x 



X e ^ 



F(0) 



3(1 - SukfHj 



(// - A)l]fc„ ^_2, _ (2/x - 3(1 - d^)X)nmo 



li-X-1 



2{ii - 1) - 3(1 - d2)A 



(27) 



(28) 



Consequently, the relation x = gives the evolution 



V2 
Ho 



4(/x-A)l],„ 2. _ (2/x - 3(1 - (j^)A + 2d^)nmo sa-d^) 
3(//-A-l) 2(/x - 1) - 3(1 - c|2)A 



b(6X — 2a + 2) -2(^-1)3: 
+ 3A " ' 



2(/i-A)rifc 2x /o,. Q^i J2^ 



(//-A-1) 



- (2/i-3(l -d^)A 



+ 6d2)(2(// - 1) - 3(1 - c/2)A)-iQ„„e-=^(i-'^')^ + -^(4A - 2// + 2) 

A 



-2{pt-l)x 



X e ^ 



fco 



-2a; 



ili-X-l) 



2n^,e-'('-''>{2{i,-l) 



3{l-(f)X)-^ + be^^ 



(29) 
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Figure 8: Plot of versus x for K-essence model. 
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Figure 9: Plot of x versus x for K-essence model. 
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Figure 10: Plot of V versus for K-essence model. 
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From here we can plot cj) and x versus x with previous assumptions. We see 
from Figure 8 that the scalar field (j) increases. Figure 9 shows that x almost 
lies in the required interval (|, |) from early epoch to the later time. The 
EoS parameter fl26l) indicates that the accelerated universe can be obtained 
for this interval. For x < |, it gives phantom DE era which corresponds 
to late time attractor [SH]. Figure 10 shows that V{(f)) increases slowly but 
attains very large negative value with the increase of scalar field 0. 



3.4 New Holographic Dilaton Field 



The Lagrangian of dilaton field can be expressed in terms of pressure of scalar 
field as [20] 

Pd = -X + &ie'^V, (30) 

where hi and 62 are taken as positive constants. This Lagrangian (pressure) 
produces the following energy density 



Pd 



-X 



(31) 



The EoS for dilaton DE is 

_ -1 + hie^-"t>x 
'^'^ ~ -1 + 36i e^'^^x 

By setting = Pa and pd = Pa, we obtain 



(32) 



e'"^X = T- 

Oi 



3(/x-A-l) 



{2p-'i{l-(P)X + 2(P)n 



2{p - I) -?>{!- 



rm_^-J,{l-d?)x 



b(6X — 2n + 2) -2(t,-i}x 
+ 3A 



2(/i- A)l]fco 



-2x 



Q. e 

^ '"mo ^ 



-3(l-d2)a 



(2// - 3(1 - (i2)A + 6^2) b.^^ ^ ^. -2(.-i). 
" 2(,-l)-3(l-dV ^A(^'"'"^'^^ ' 



(33) 



The EoS parameter ( !32l) gives the bound of e^'^'^x which is (y, y) in 
order to obtain the accelerated universe. Its graph versus x with a = 0.05 is 
shown in Figure 11 which indicates that this almost lies in the same interval 
and shows consistency. Also, the solution of (133|) follows 
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Figure 11: Plot of e^^^x versus x for dilaton field. 
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1^0 



4(/j, - \)Vtko 2x _ Q 3(l-d2)x 

3(//-A-l) 



i2iL-m- (P)X + 2(P) 6(6A-2/x + 2) ^2(M-i). 

X -^TT^T rr- 3TTT^ H ^ rr-^- ^ 



X 



2(2(/i - 1) - 3(1 - t;2)A) 3A 
2(/i-A)r],, _ (2//-3(l-(i^)A + 6rf^)a^, _3(,_,.), 
(^-A-1) 2(;u - 1) - 3(1 - d2)A 



+ -(4A-2/i + 2)fee^^ 
A 



1] 



-2z 



(/x-A-1) 

X (2(^ - 1) - 3(1 - (f)X)-'e-'^^-'''^^ + 6e"^ 



2fi, 



dx 



(34) 



Its graph is shown in Figure 12 which exhibits direct proportionality with 
respect to x. This type of behavior gives scaling solutions for dilaton model 
as proved in [20]. 



4 Generalized Second Law of Thermodynam- 
ics 

In this section, we check the validity of GSLT for interacting NHDE in the 
non-fiat universe. Bekenstein's [H] provided a relationship about the en- 
tropy of black hole horizon and horizon area which plays a crucial role to 
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Figure 12: Plot of versus x for dilaton field. 



discuss the GSLT. This law states that the sum of black hole entropy and 
the background entropy must increase with the passage of time. The first 
law of thermodynamics yields 

TdS = pdV + dE, (35) 

where T, S, p, V and E denote temperature, entropy, pressure, volume and 
internal energy of the system, respectively. Differentiating it with respect to 
time, we obtain 

• paV + Ex ■ PmV + E^ 

- J, , - J. , [OO) 

for the NHDE and DM, respectively. The volume, temperature and entropy 
of horizon L in non-flat universe are defined as 



47^ 5^ = 27r2L2. (37) 

3 ' 27rL' ^ ' 



The internal energies for NHDE and DM are 



£, = 1^, = (38) 



Using Eqs.(I7]), (|9]) and (l36!) - (!38l) . we have the final expression of the GSLT 

/I 2 2 / 

SLai = + ^a)Pa + 3H',n^,e-'^'-''y^]{p'j, + 2pA) - 

17 



1.0 

0.5 
_ 0.0 
-0.5 
-1.0 
-1.5 

0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

X 

Figure 13: Plot of S^^^^i versus x for NHDE. 

where S'^^^^i denotes the derivative of the sum of entropies due to DM, NHDE 
and horizon entropy. We plot it against x with the same assumptions as given 
earlier. Figure 13 shows that the GSLT is valid for early epoch but fails for 
the later time in this scenario. 

5 Concluding Remarks 

It was shown through experimental data that our universe is not totally 
flat but it contains small positive curvature [26]- [M]. This motivated us 
to perform the versatile study of the interacting NHDE in the context of 
non-fiat universe. In this paper, we have discussed four main features by 
choosing different values of interacting parameter i.e., d"^ = 0, 0.04, 0.08 which 
are summarized as follows. Firstly, we obtain the evolution equation which 
starts from vacuum DE region {u\ = — 1) and goes towards quintessence 
region {u\ > —1) as shown in Figure 1. This type of behavior is consistent 
with the present observations. Secondly, we have discussed the instability 
of this model against perturbation as shown in Figure 2. It is found that it 
remains unstable forever like HDE with IR cutoff as a future event horizon 
[17] . agegraphic [18] and ghost QCD [19] DE models but unlike Chaplygin 
gas [16] DE model. 

Thirdly, we have evaluated the interacting NHDE versions of quintessence, 
tachyon, K-essence and dilaton scalar field DE models. We have also explored 
the behavior of potentials and dynamics of the scalar field corresponding to 
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these models by using the rehable values of the constant parameters. These 
are shown in Figures 3-12. It is seen that the scalar field <^ becomes more 
steeper with the increase of the interacting parameter cP for all scalar field 
models. Also, the scalar potential V{(j)) decreases more rapidly in case of 
quintessence for increasing the interacting parameter rf^ as compared to non- 
interacting scenario (Figure 5). However, for K-essence model, it increases 
but achieves a maximum value lower than that of non-interacting case as 
shown in Figure 10. It was argued [2T| H3l HI] that the scalar field models 
favor the DE phenomenon for EoS parameter lying in the interval (—1,0). 
However, in our case, the EoS of NHDE remains in this region and shows 
consistency. 

We would like to mention here that our solutions coincide with the attrac- 
tor solutions (for quintessence [21] and K-essence [S^ DE models) and scaling 
solutions (for tachyon [38] and dilaton [20] DE models). It is remarked that 
our expressions of EoS parameter and scalar field models can be reduced to 
the results of [22] (with vanishing of DM density and A; = 0) and of [23] (for 
the vanishing of DM density only). Finally, we have checked the validity of 
the GSLT in this scenario. It is found that initially it is valid but fails for 
the later time. 
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